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The measurement of the flavor composition of the neutrino fluxes from astrophysical sources has 
been proposed as a method to study not only the nature of their emission mechanisms, but also 
the neutrino fundamental properties. It is however problematic to reconcile these two goals, since a 
sufficiently accurate understanding of the neutrino fluxes at the source is needed to extract informa- 
tion about the physics of neutrino propagation. In this work we discuss critically the expectations 
for the flavor composition and energy spectrum from different types of astrophysical sources, and 
comment on the theoretical uncertainties connected to our limited knowledge of their structure. 
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I. INTRODUCTION 

There is the expectation that in the near future we will see the opening of the new field of observational high 
energy neutrino astronomy 

0, 0, II ■ Theoretically there are very robust reasons to expect the existence of high energy 
(E v > 10 12 eV) neutrino sources. The strongest motivation is the observation of a cosmic ray flux, mostly composed of 
protons and fully ionized nuclei that extends in energy up to E ~ 10 20 eV. These hadronic particles can interact inside 
or near their acceleration sites or during their propagation in interstellar or intergalactic space. These interactions 
produce weakly decaying particles (such as tt^ and kaons) that generate neutrinos. These "astrophysical neutrinos" 
are intimately connected with the high energy photons that are created in the decay of 7r° and r\ particles produced 
in the same hadronic primaries interactions, or in the radiation processes of relativistic electrons and positrons co- 
accelerated in the sources. A rich variety of high energy gamma ray sources has been observed with detectors on 
satellites [3| and ground-based 0] Cherenkov telescopes, suggesting several possible neutrino sources. 

Observations with neutrinos have the potential to give us unique information about their astrophysical sources, 
and hopefully could also result in the discovery of new classes of sources. On the other hand the possibility to use 
these observations to obtain information about the fundamental properties of the neutrinos has been widely discussed 
0, d, H, [l(| [HI [H, EH- Astrophysical neutrinos travel pathlengths of order 10 Kpc (~ 3 x 10 22 cm) for galactic 
sources, and as large as ~ 1 Gpc (~ 3 x 10 27 cm) for extragalactic sources. These remarkably long baselines allow the 
study of phenomena such as v flavor transitions or v decay in a range of parameters that is unaccessible with other 
methods. Only observations with SuperNova neutrinos, very likely only observable from galactic sources, but having 
much smaller energy {E v ~ 10 MeV) could provide larger L/E v . 

In order to extract information about the neutrino fundamental properties from the observations one needs to have 
sufficiently good understanding of the properties of the neutrino emission at the source. The common assumption in 
several studies of this type is that the properties of the v emission, and in particular the flavor composition at the 
source can be robustly predicted. There is of course a well known historical precedent for the successful use of this 
concept in the discovery of neutrino oscillations with atmospheric neutrinos. In fact, the first hints of the existence 
of the now solidly established «-> v T transitions were obtained with the Kamiokande [f| and 1MB Q detectors as 
the measurement of a fi/e ratio for contained events smaller than the expectations. It is natural to try to make use 
of new, very distant neutrino sources (when they will be discovered) to perform additional studies of the properties 
of neutrino propagation. 

Several "exotic" processes, beyond standard flavor oscillations, could reveal themselves only in the propagation of 
neutrinos over very long distances. For example, it has been suggested that some neutrinos could decay into a lighter 
neutrino and a majoron [lij], if the lifetime is sufficiently long this phenomenon could be only detectable for neutrinos 
propagating over astronomical distances (l5l fl6| . A second interesting possibility is that neutrinos are pseudo-Dirac 
states [17] where each generation is actually composed of two maximally-mixed Majorana neutrinos separated by a 
tiny mass difference. If the pseudo-Dirac splittings are sufficiently small, the phenomenology of oscillations on short 
baselines remains unchanged, however when E v j L becomes comparable or smaller than the pseudo-Dirac splittings 
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new transitions become possible, and can in principle be detectable with astrophysical neutrinos [18j. More in general, 
oscillations into sterile states that are quasi degenerate to the active neutrinos can in principle be investigated down to 
very small squared mass splittings [HI . Several other mechanisms such as quantum decoherence [2(| , violations of the 
equivalence principle [U, [22j , neutrinos with varying masses [23|, [24[ could leave their signature on the propagation 
of astrophysical neutrinos. 

A well known "naive" argument states that since the dominant source of astrophysical neutrinos is the decay of 
charged pions, and each n , after chain decays such as: ir + — > /i + z/ M — > (V^ e + v e ) v^, generates two muon neutrinos 
and one electron neutrino, the flavor ratio at the source is R^ e = (y^ +V fl )/(i> e + V e ) ~ 2. This naive argument 
has some pedagogical value, but must be considered only as a first approximation. In fact, even if charged pions are 
the only source of neutrinos, i? Me is only equal to 2 after integration over all v energies, because the three neutrinos 
generated in a charged pion decay have different energy spectra, and therefore in general the flavor ratio will depend 
on the spectral shape of the v signal, and will vary with the neutrino energy. An effect that can be of large importance 
is the presence of sufficiently efficient energy loss mechanisms in the source. In this situation the particles that are 
most affected are the muons. The presence of additional neutrino sources (such as kaons) can also be important for 
the flavor ratio. In this work we will re-examine critically the uncertainties in the predictions of the spectra and flavor 
composition of astrophysical neutrinos and discuss the implications for the extraction of information on the neutrino 
properties. 

This paper is organized as follows: in the next section we discuss how the observable neutrino flavor ratios carry 
information at the same time about the flavor composition at the source and about the flavor transition probabilities. 
In section 3 we outline the general structure of the calculation of the neutrino signal from an astrophysical source. The 
different steps of these calculations are considered in more details in sections 4-7. In section 8, as an illustration, we 
describe an explicit calculation of the neutrino emission from the fireballs of Gamma Ray Bursts following the model 
of Waxman and Bahcall [25] . Section 9 contains a short discussion of the problem of the experimental determination 
of the flavor composition and energy spectra of a neutrino signal. The last section presents our conclusions. 

II. FLAVOR RATIOS OF ASTROPHYSICAL NEUTRINOS 

The observable fluxes of astrophysical neutrinos from a source at distance L will be linear combinations of the fluxes 
at the source. Leaving implicit the energy and distance dependences one can write: 

4>v a = ^ p u^v a 4>° Up ■ (!) 

The transition probabilities P„__n, a have certainly a non trivial structure because of the known existence of "standard" 
flavor oscillations, but might depend on additional "new physics" contributions, such as neutrino decay, that become 
significant only for very long pathlengths. In the standard scenario the neutrino number is conserved and therefore 
^2pPa->p = 1- More in general, in the presence of a non negligible decay probability or of transitions to additional 
sterile states, the sum can be less than unity. 

The flavor, energy and pathlength dependences of the standard oscillation probabilities are well known, and are 
determined by two squared mass differences (AmJ 2 ~ 8.0 x 10~ 5 eV 2 , |Am| 3 | ~ 2.5 x 10~ 3 eV 2 ) and the neutrino 
mixing matrix. The oscillation lengths (Ay = 47ri? ! ,/|Am 2 J) of the standard flavor transitions are short with respect 
to the typical astrophysical distances and, in most cases, it is a good approximation to consider only the probability 
averaged over distance. In this case the probability becomes independent from E v and L and takes the form: 

3 

where U is the unitary mixing matrix that relates the neutrino flavor and mass eigenstates. This matrix can be 
written in terms of three mixing angles and one CP violating phase. From a global fit to all existing neutrino data 
[26| one can extract the best fit values: 0\2 — 34°, 6*23 — 45°, 6*13 ~ 0, and 99% C.L. intervals for the mixing angles: 
6»i2 € [30°, 38°], 6» 23 S [36°, 54°], 6 13 < 10°; the phase S remains completely undetermined. 

Assuming to know sufficiently well the properties of the source, the observed flavor ratios can give information on 
the flavor transition probabilities. This can be used to help in the determination of the "standard" parameters in 
the flavor oscillations, or more ambitiously to investigate the possible existence of additional phenomena in neutrino 
propagation. 

As an illustration in fig. [1] and [2] we show the expectations for the two independent flavor ratios R e u = (ve/^^obs 
and R^r = (^,/^ r )obs, calculated for different assumptions for the source emission, in the presence of simple standard 
oscillations, or also including neutrino decay (with two different assumptions) [161 ] . The distributions of the flavor 
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ratio for a given model are determined only by the present uncertainties on the neutrino mixing parameters. A 
source model is defined by the relative intensity of the emission for the three neutrino flavors. For our illustration 
we have considered three source models. The first model is the emission of form: [i/ e , is u , v T \ SOU rce = [1,1.86,0]. This 
corresponds (as we will discuss more extensively in the following) to pion dominated emission from a thin source (no 
significant energy loss for secondaries) with a power law energy spectrum of slope 2. Several classes of v sources of 
this type (such as Supernova Remnants and Gamma Ray Bursts) have been predicted. Note that the flavor relative 
abundances in this model are close but not identical to the naive expectation [1,2,0]. The deviations are a simple 
consequence of the shape of the energy spectra of the neutrinos in pion chain decay. The observable R efl ratio for 
this model takes value in the interval ~ (0.7, 1.2) with a most likely value close to unity, while the other independent 
ratio R UT has a narrow distribution sharply peaked at the value 1, with a tail extending up to ~ 1.3. 

Figures Q] and [5] also show the predicted flavor ratios for emission with (summing over v and 77) flavor abundances 
[0,1,0], and [1,0,0], that is pure muon or electron neutrino emission. The motivation for using these two models 
is that, to a very good approximation, they can be considered the two extreme models for the emission from a 
"standard" (not involving new physics) source, since a significant production of v T in an astrophysical environment 
is extraordinary unlikely. Sources emitting pure fluxes of v a or v e are in principle possible, and have been in fact 
advocated in the literature. They correspond to a source dominated by pions where muons lose all their energy before 
decay, or to a pure neutron source, that generates a v e flux. A pure neutron (T> e ) neutrino source is in fact unrealistic 
since it is natural that neutron production is associated with pion production and therefore with some v a emission. 
We will comment more extensively on these issues in the following. 

For the [0,1,0] ([1,0,0]) model the most likely value for the observable flavor ratio is R efl ~ 0.55 (R ef i — 2.7). In both 
cases R aT is centered at the value unity, this is a simple consequence of the large (approximately maximal) mixing 
between v a and v T . In the case of a [1,0,0] source the R^ T distribution is noticeably broader extending to the interval 
RfXT G (0.5, 1.5). 

Neutrino decay with a very long lifetime is one of the most interesting possibilities that can be investigated with 
astrophysical neutrinos. Assuming [13, EE] that only the lowest mass eigenstate is stable, and that the distance of the 
source is much longer than the decay lengths, the observable flavor ratios depend uniquely on the flavor composition 
of the lightest eigenstate. Depending on the sign of Am§3 this is v\ (direct mass hierarchy) or i/ 3 (inverse hierarchy). 
The flavor ratios for these two models are also shown in fig. [T] and [2] as dashed lines. The R eiL ratio is particularly 
interesting. In the case of direct mass hierarchy, the most probable value of the e/j, ratio is large (R eu ~ 4.2), the 
distribution is also remarkably wide, allowing in principle to constrain the mixing parameters [161 ]. For the inverse 
hierarchy, when the stable neutrino is the v% that has a small (or perhaps vanishing) overlap with v e , the observable 
R eu is close to zero. 

It is remarkable that the predictions for the e/i flavor ratio in neutrino decay are in both (direct and inverse mass 
hierarchy) models significantly different from the standard model predictions for nearly all assumptions about the 
flavor composition of the emission. Only a pure v e emission has some overlap with the decay model with stable v\. 
The study of the e/i flavor ratio of astrophysical neutrinos can therefore give evidence in favor or against the existence 
of v decay. 

The determination of the neutrino mixing parameters in the standard oscillation scenario is more problematic. For 
the most sensitive e/i ratio, the range of possible values due to the present uncertainties on the mixing parameters is 
in fact smaller than the variations that result from different flavor abundances at the source. Therefore the possibility 
to obtain interesting bounds on the mixing parameters depends crucially on having a sufficiently precise knowledge 
of the source. 

To illustrate this problem, we can write an approximate expression for the flavor ratio R e ^ in terms of the relevant 
mixing matrix parameters expanding in first order around the best fit values 

[1,1.86,0] 
obs 



1.026 - 0.0218 80° 3 + 0.0087 6° 13 cos 8 - 0.0014 56° 12 + 0.370 8 — (3) 



\ [04,0] / 

— J ~ 0.547 -0.0296 5023 + 0.0119 13 cos 8 + 0.0197 86^ + 1-152 8 (— ) (4) 



s [i»o,o] 

v, \ 1 _„„ , „ „ / v. 



2.65 + 0.093 (56*23 - °- 037 9 i3 cos s ~ °- 131 6d i2 ~ 3-846 <M — I (5) 

V V /obs V v e J o 

In these equations 80\ 2 and 89^ are the deviations in degrees of the mixing angles from their best fit value (6*12 = 34°, 
#23 = 45°), # 13 is the value (always in degrees) of the third angle and 8 is the CP violating phase. The superscript 
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labels of the flavor ratio indicate the source model. The last term in equations ([3]), |@J and ([5]) gives an estimate of 
the variation in the observable ratio due to the uncertainty in the flavor ratio at the source obtained expanding in first 
order in 8{v a /vfj)o around the flavor abundances of the model considered. One can see that uncertainties in flavor 
abundances at the source can be so large to make impossible a meaningful constraint on the mixing parameters. 

In the remaining sections of this paper we will discuss how the properties of the astrophysical sources determine 
the flavor abundances of their emission, and estimate the associated theoretical uncertainties. 

III. PRODUCTION OF ASTROPHYSICAL NEUTRINOS 

Astrophysical neutrinos are generated when a population of relativistic hadrons (protons or nuclei) interacts with 
some target material (gas or a radiation field) inside, near, or outside the acceleration site. These interactions produce 
weakly decaying secondary particles whose decay generate neutrinos either directly, or indirectly with the subsequent 
decay of muons. The description of these processes therefore requires the following elements: 

1. The description of the energy spectrum and composition of the primary particles. 

2. The definition of the target material with which the primary particles interact. 

3. The modeling of the properties of particle production in hadronic interactions. 

4. The description of the properties of the medium where the interactions are taking place, to determine the 
relevant mechanisms for energy loss. 

5. The (well known) properties of weak decays. 

In our discussion we will describe the source as one stationary homogeneous region. This is an important limitation, 
because in general we can expect that the sources will have non trivial space structures, with the neutrinos emerging 
from different regions that can produce different spectral shapes and flavor contents; most high energy radiation sources 
are also expected to show important time variabilities. Because of the limited statistics and angular resolution, the 
observations of astrophysical neutrinos will necessarily integrate over the entire volume of the source (that in most 
cases will appear as point-like) and average over most or all time variations. Describing this situation requires the 
averaging over different source regions and conditions. 

Some of the most interesting postulated sources of neutrinos, in particular the jets of Active Galactic Nuclei (AGN) 
and Gamma Ray Bursts (GRB), are expected to be associated with astrophysical jets with ultra-relativistic bulk 
motion. These relativistically moving sources are best described in the rest frame of the jet, where the emitted 
radiation is approximately isotropic, and the electromagnetic fields can be described as a purely magnetic field. The 
observable fluxes can then be obtained by an appropriate Lorentz boost. 

The v emission from a given source can be described by the emissivities Q Va (E v ,£l) (in units (s sr GeV) -1 ) that 
give the number of neutrinos of type a emitted with energy E v in the direction 12. The v fluxes at the Earth are 
obtained from the neutrino emissivities as: 

<M^) = (1 4 ^2 )2 E p ^^) Qu,[Eu{l + Z)M (6) 

where is the source luminosity distance, z is its redshift, and P Up ^, Va are the flavor transition probabilities. The 
observable fluxes depend only on the emission in the direction 12® that corresponds to the line of sight from the source 
to the Earth. In the following we will not consider the angular dependence of the source emission, since this is not 
observable for any individual source. 

The calculation of the neutrino emissivities Q Va (E v ), starts from a decription of the energy spectrum and com- 
position of the primary particles. In the following we will assume that the primary particles are protons and are 
described by the energy spectrum N p (E p ) (in units GeV -1 ). The extension to the case where there is a significant 
contribution from nuclei is straightforward, and to a very good approximation can be calculated introducing effective 
proton and neutron spectra (N p (E p ) and N n (E n )) that take also into account the bound nucleons. The description 
of the environment of the source must include the density and properties of the target material and of all other fields 
that can be a source of energy loss. The interactions of the primary particles result in the emission of secondary 
particles , n~ , K + , . . .) that are sources of neutrinos. These particles must be propagated in the source medium 
until they decay. All muons generated by these decays must also be propagated and their decay studied. 
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The scheme of the calculation is therefore the following. The rate of production of secondary particle of type a can 
be calculated with the convolution: 

/d 
dE p N p (E p ) K p (E p ) -J^(E a -E p ) (7) 

where K p (E p ) is the interaction probability per unit time of a proton of energy E p , and dn p —> a /dE a (E a ; E p ) is the 
number of particles of type a in the final state in the energy interval (E a , E a + dE a ). 

The energy distribution of particle a at decay Q^ ec (E) will in general differ from the distribution at production and 
can be obtained with the integration: 

f°° dn a 

Qa C (E) = J dE t Q a (Ei) J^.{E;Ei). (8) 

where dp^ ec /dE(E; Ei) is the probability density that particle a created with energy Ei will decay with energy E. 
The neutrinos created in the direct decay of particle a can be calculated as: 

Q a - rUa {E„)= dE a Qf c (E a ) J£pL(E v ;E a ) (9) 

folding the parent energy distribution (at decay) with the appropriate decay spectra dn a ^ w /dE v (E v ;E a ). 

The muon production rate can be calculated with an expression that is the analogous of ((9]) substituting the 
appropriate decay distributions: 

/•OC 1 

Qm(*W=E / dE a Qf c {E a )^-^{E,;E a ). (10) 

The sum runs over all weakly decaying hadrons. 

In out calculation we explicitely included charged pions, kaons and neutrons, neglecting the much smaller contri- 
butions of heavier particles. An important complication is that the spectra of the particles generated in muon decay 
depend on the muon helicity (see discussion in section Hv|) . It is therefore convenient to calculate separately the pro- 
duction of muons with different helicity: {//J, /z^, /z^ and /z^} (where L and R indicate the helicity). It is expected 
that the helicity of the muons is to a good approximation conserved even in the presence of a strong magnetic field 
and of significant energy losses, and therefore the energy distributions of the muon of different helicities at decay can 
be obtained with a convolution similar to (151): 



r°° dr^ 
Qj(£) = J e dEi Q /J ±(E i ) -^^(E;Ei 



(11) 



The last step is to compute the neutrinos generated by muon decay: 



r°° dn n ±^„ 
Q^ V JE V ) = J e dE, Q d ^(E,) ^ ° (g„;ff M ). (12) 

The neutrino flux is obtained summing over all possible sources. 

In the following sections we will give more details about the neutrino production. Our discussion starts from the 
neutrinos and goes backward along the chain of processes that lead to their production. The next section contains a 
discussion of the decay energy distributions; after discussing the consequences for exact power law spectra in sectionlVl 
in section [VI] we will discuss possible mechanisms for energy loss and the calculation of the decay energy distribution; 
in section rVlII we will discuss the interaction probability K p (E p ) and the properties of particle production in hadronic 
interactions. 



IV. ENERGY SPECTRA IN WEAK DECAYS 



For an accurate prediction of the neutrino fluxes it is necessary to include a precise description of the energy spectra 
produced in weak decays. The energy distribution of particle b in the decay a — ► b + . . . in the frames where the 
parent particle a is ultrarelativistic, takes the scaling form: 

''"- b (E b] E a ) = ±F a ^(^). (13) 



dE v E a \E a 
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The functions F a ^b(x) for the relevant weak decays are calculable using the measured branching fractions in the 
possible final states, and their matrix elements. 

For the two body decay 7r + — > that accounts for approximately 100% of charged pion decay, the scaling 

distributions are fully determined by elementary kinematics: 

F 7r+ ^(x) = —!—0(l-ri T -x) (14) 

F v+ ^+ (x) = —— G(x - rv) (15) 

where 0(x) is the step function and r v — (m^/m^) 2 . It is also necessary to take into account the spin state of the 
final state muons, that can be described by the helicity h = Pr — P^, where Pr.l is the probability that the muon 
has spin parallel (anti-parallel) to its momentum. The helicity, in a frame where the parent pion is ultrarelativistic, 
is a function of the fractional energy x [13] : 

(x) = ~ TTT' (16) 

Because of CP invariance: 

K-^^-ix) = -K+^+ix) (17) 

The helicity h 7r +_ >fi +(x) takes the value —1 for x — > 1 (forward muon emission) and +1 for x — > r n (backward 
muon emission) reflecting the fact that the fi + is created as a left-handed particle in the pion rest frame in order to 
compensate the angular monentum of a left-handed neutrino. 

A useful method to take into account the effect of muon polarization in the presence of energy loss for the muons is 
to consider separately the production of left-handed and right-handed muons, that have two well determined decay 
spectra. Combining equations (fT5]) and (HHJ) one obtains the distributions: 



F 



\ + ^ t {x) = F„-^ +R {x) = {l _ r J )2x eCar-r,) (19) 

Because of CP invariance one has: F + + (x) — F _ - (x). 

The above discussion is also valid for the decay mode K + — > (and charge conjugate), with the simple 

replacement m, — > vtik- 

The scaling functions that describe the decay of a muon of helicity h are: 

F^^Jx;h)=^-3x 2 + ^+h (-I + 3z 2 -^-) (20) 

Fp+^ve {x; h) = (2 - 6 x 2 + 4 a; 3 ) + h (2 - 12 x + 18 x 2 - 8 x 3 ) (21) 

Because of CP invariance the spectra for the charged conjugate decays are obtained with the replacement h — > —h. 
The spectra for unpolarized muons can be obtained setting h = 0. 

If the muon energy loss before decay is negligible, it is straightforward to convolute the previous expressions to 
obtain the v spectra after a chain decay (ir — > fi — > v) or (K — > fi — > v). In the more general case the energy loss of 
the muon must be taken into account. 

The energy spectra of neutrinos and muons emitted in the three body decays of kaons (K e $° and K^°) can be 
written in terms of vector form factors that have been experimentally determined [28| . 

Neutrons are also a source of T> e . The energy spectrum of these neutrinos is also well known. The v e carry a very 
small fraction of the parent neutron energy and are negligible in most (but not all) cases. 

Several works on the fluxes of astrophysical neutrinos make simplifying assumptions about the energy distributions 
of the neutrinos produced in weak decays. For example the calculation of Kashti and Waxman [2^| approximates the 
decay spectra of charged pions and muons into muons and neutrinos as: 



Fn+^Vu, ( x ) = S 



1 

X "4 



(22) 
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(23) 



F, 



fl+ — >V, 



(x)=S 



1 

X ~3 



(24) 



These approximations, together with the assumption that muon energy loss before decay is negligible, imply that a v 
source dominated by pions has a flavor ratio [y^ + V^)/(v e +V e ) = 2. With the use of the correct expressions (|18|) . 
(TT9)l . (|2"D|) and (|2"Tj) for the decay spectra one obtains the result that even for a pure pion source the flavor ratio is not 
exactly two, and is in general a function of the v energy determined by the shape of the neutrino spectrum. 



V. POWER LAW SPECTRA 



The situation where the neutrino energy spectrum is a power law of form <^>„ cx E~ a is phenomenologically very 
important, and is simple to discuss. A power law v spectrum implies that the spectrum of the parent particles (pions 
and kaons) is also a power law with the same slope (a Wt K — a u)- Such a spectrum arises when the interacting primary 
particles have a power law spectrum and the target is either composed of normal matter at rest, or is a photon field 
with an energy distribution that has again a power law form. In the first case the slope of the secondary particles 
spectrum is equal to the one for the primary particles (a Vt K — ot p ); in the second case the slope of the secondaries 
is a^.K ~ a p — f3 + 1 where a p the slope of the primary particles and [3 the slope of the target photon distribution 
(n 7 oc e^ 13 ). 

In general the spectrum of neutrinos produced in the decay of particle type a can be obtained convoluting the 
energy spectrum of the parent particles with the appropriate weak decay spectrum (see equation @). If the parent 
particle spectrum has a power law form, because of the scaling form (|13p of the decay spectra the resulting neutrino 
energy distribution can be written as: 

f°° dn 
Q V {E V ) = / dE a Q a (E a ) -£pL(E v ;E a ) 
Je v dE v 

= C a E- a f dx x^ 1 F a ^(x) = C a Z a ^{a) E- a . (25) 
Jo 

In other words the neutrino spectrum is a power law with the same slope of the parent particles, with a proportionality 
factor that is commonly called the "Z-factor" , and is the momentum of order a — 1 of the decay spectrum: 

Z a ^ 3 {a) = I dx Z" 1 F a ^ o {x). (26) 
Jo 

Z a ^ v . (1) is the average multiplicity of Vj in the final state, and Z a ^ v . (2) is the fraction of the parent particle energy 
carried away by neutrinos of type j. 

In the case of tt^ chain decay the three relevant Z-factors are: 

Z w+ ^(a)= [ - (27) 



4 (3-2rv-(3 + a ) rg + (2 + a ) rj+ a ) 
a 2 (2 + a) (3 + a) (1 - r^f 

M a 2 (1 + a) (2 + a) (3 + a) (1 - r x ) 2 

with r w = (m^/m^) 2 . The Z-factors for charged conjugate modes are identical because of CP invariance. The 
Z-factors given above are calculated neglecting the energy loss of muons before decay, and assuming that the muon 
helicity is exactly conserved. This assumption remains valid, to a very good approximation, also in the presence of a 
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magnetic field, because the bending of the momentum and the spin precession exactly cancel for a particle of electric 
charge ±1 and magnetic moment of one Bohr magneton. Neglecting the effects of the muon polarization leads to 
an overestimate (underestimate) of the Z-factor for the ir + — > /i + — > F M (tt + — ► n + — > v e ) channel. For a = 2 the 
Z-factors for pion decay take the values: 

S7 7 7 \ f LlIl 3 + 4r ^ 2 + M 

i'"7r+— >v )1 ) •^7T+— >^i+— >i/ M ) ^tv+— *fi+— >i/ e /a=2 > 2 ' 20 ' 10 I 

(30) 

~ {0.2135, 0.2646, 0.2573} , 

i.e. the three neutrinos carry approximately one quarter of the charged pion energy. This happens because in the 
first decay the muon carries away a large fraction ((1 + r 7r )/2 ~ 0.787) of the pion energy. 

Because of the different shapes of the energy distributions of the three neutrinos emitted in a pion decay, the flavor 
ratio 

-K/xe = 7 — r = r (31) 

(y e + v e ) Ar+-»/i+-n/ e 

for a power law spectrum of parent pions is a function of its slope. The ratio is shown in fig. [3J For a = 1 one has 
i? Me = 2, since in this case the Z factors correspond to the neutrino multiplicities. With increasing a the flavor ratio 
decreases monotonically, this reflects the fact that the muon neutrinos produced in the direct pion decay are softer 
than the neutrinos of the muon decay. For the value a = 2 one has i? Aie ~ 1.858. Figure [3] also shows the predictions 
obtained neglecting the effects of muon polarization to illustrate the importance of their inclusion. 

The second most important source of neutrinos is the decay of kaons. Charged kaons can produce neutrinos with 
three different decay channels: the two body mode K + — > [i^v^ (with branching ratio 0.6343), and the three body 
modes K% (K+ -> n°e+v e ) and K± 3 (K+ -> 7r> + i/ M ) that have branching ratios Bf 3 ~ 0.0487 and B± 3 ~ 0.0327 
[28j ] . The Kl can also produce neutrinos in the K° 3 decay mode (Kl — * 7r =F e ± ^ e (77 e )) with combined branching ratio 
B® 3 ~ 0.3881, and K° 3 mode (Kl — * 7r =F ^ ± ^ AI (77 M )) with combined branching ratio B® 3 ~ 0.2719. The i? Me ratios for 
neutrinos produced by the chain decay of charged or neutral kaons is shown in fig. [41 In the case of charged kaons 
the 2-body decay mode is the dominant one, the Z-factors for this mode have the same form as in pion decay, with 
the replacement r n — » rjc, however for charged kaons the flavor ratio grows when a increases, reflecting the fact that 
the muon neutrinos produced in the direct decay carry nearly half ((1 — r^-)/2 ~ 0.48) of the parent particle energy 
while the neutrinos produced in the second stage (from muon decay) take on average only 18% (u^) and 16% (y e ) of 
the kaon energy. The inclusion of the three body decay modes (K^ 3 and K^ 3 ) that have smaller branching fractions 
reduces the flavor ratio, but it is a small correction. 

In the case of the decay of Kl, the three body decays are the only source of neutrinos. The K° 3 mode has a 
larger branching fraction than the K° 3 channel, because of the larger phase space available. When a = 1, when the 
flavor ratio reflects the neutrino multiplicities in the final state, the flavor ratio is i? Me = 2B° 3 /(B° 3 + B°^ 3 ) ~ 0.82. 
Increasing a the ratio decreases, because the v e from the K° 3 decay have the hardest spectrum, and their contribution 
is enhanced. 

Note that for power laws with slope ~ 2 the net effect of including kaon decay as a v source is a only a small 
(positive) correction, because of a cancellation between the contributions of charged kaons (that increase the relative 
importance of muon neutrinos) and neutral kaons (that do the opposite) . 

VI. ENERGY LOSS MECHANISMS 

In the most general case, the unstable particles that create neutrinos can lose a significant amount of energy during 
the time that elapses from the moment of their creation to the moment of their decay. This effect can have very 
important consequences for the spectrum and flavor composition of the neutrinos. 

As an illustration, muon energy losses are very important for the prediction of the atmospheric neutrino fluxes. The 
largest effect is due to the fact that the muons that reach the ground lose rapidly energy because of ionization and 
electromagnetic radiation processes, and decay practically at rest (or are captured by heavy nuclei [30l ]) producing 
only very low energy neutrinos. Pion energy losses are also very important for the calculation of the atmospheric v 
fluxes. In this case the dominant mechanism for energy loss are hadronic interactions in air. The pion decay length 
grows linearly with E„, and (for E v ^ 100 GeV) it becomes (for vertical particles) comparable to the interaction 
length. Most pions of higher energy interact before decaying. These hadronic interactions result in a multiplication 
of the pions, but the secondary particles have much lower energy and the net effect is a strong suppression of the 
neutrino flux at high energy. 
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There are several mechanisms for energy loss that could be present in the astrophysical environment where neutrinos 
are produced. In several circumstances synchrotron radiation could be the dominant source of energy loss [29l l3l| . In 
the presence of a magnetic held of value B, the energy loss of a particle of electric charge e, mass m and energy E, 
after averaging over all possible orientation of the magnetic held direction is given by the well known expression: 

dE Ae i B 2 „ 

~-r = t: — a- E (32) 
dt 9 to 4 v ' 

The energy loss is significant only for E > E syn the critical energy where the the synchrotron loss time i syn = 
— [{dE / dt) / E\ _1 is equal to the decay time idee = t E/m: 

3 m 5 / 2 _5.8xl0 18 eV / m \ ^ f r \~* 

The critical energy for synchrotron losses scales oc to 5 / 2 t -1 / 2 . It is smallest for muons, it becomes 18.4 times larger 
for charged pions, and 628 times larger for charged kaons. There is therefore an energy range where the losses are 
only significant for muons, and a second energy range where the synchroton losses are significant for both pions and 
muons but not for kaons. 

In the presence of ordinary matter in gas form the particles can also lose energy because of ionization and radiation 
processes (with radiation dominating at large energies). These energy losses take the form —dE/dt ~ p(a + bE) 
where p is the medium density and a and b are slowly varying coefficients that depend on the gas composition. For a 
medium mostly composed of protons, as it is likely in astrophysical environments, the critical energy for the energy 
losses with ordinary matter in gaseous form (where the loss time equals the decay time) is for muons: 

Matter - 6.5 x 10 18 (p-io)" 1 eV (34) 

where p-io is the matter density in units of 1CP 10 g cm~ 3 . 

In the presence of a gas of ordinary matter the dominant source of energy loss for charged pions and kaons is due 
to hadronic interactions. The critical energy at which the hadronic interaction time t[ nt = (A) / (Na p Pint) equals the 
decay time r E/m can be estimated for pions as: 

^matter * 7.4 X 10 16 (p^ 1 eV . (35) 

This critical energy scales oc m/r and for charged kaons (assuming the same hadronic cross section) is approximately 
7.4 times larger. 

The effects of ordinary matter on the secondary particles can be safely neglected for all energies (including those 
above -E ma tter) and all particle types if the density is below the critical value: 

p crit ~ 2.0 x 1(T 15 B Gauss g cm' 3 . (36) 

For p < p cr jt the energy losses on matter are either negligibly small, or dominated by the synchrotron losses. 

Figures O and [5] illustrate the above discussion showing the characteristic times idec (decay) , t syn (synchrotron 
losses) and i ma tter for muons and pions. 

Gamma Ray Bursts have been proposed by Waxman and Bahcall [25j as a neutrino source having a potentially 
very strong magnetic field. In the WB model, the strength of the magnetic field is estimated by energy equipartition, 
assuming that it corresponds to an energy density of the same order of magnitude as the energy density carried by 
the fireball photons. It is natural to ask the question if this radiation field can be an important source of energy 
loss for secondary particles produced in it. The calculation of the time tjc for Inverse Compton losses of charged 
particles traveling in the fireball radiation field requires to take into account the energy spectrum of the target photons 
(assumed to be isotropic in the GRB jet frame). In fact the Compton scatterings with sufficiently high energy photons 
e > m 2 / E happen in the Klein-Nishina regime and are inefficient as a source of energy loss. For the predicted GRB 
photon spectrum [25j |. (given also in equation (|63p ) that falls at high energy oc e -2 , the energy loss for Inverse Compton 
grows at high energy as — (dE/dt)ic oc E p 1 /m 2 , (where p 7 is the target photon energy density), and the loss time 
tie becomes constant. The results of a detailed integration are also shown in fig. [5]and[HJ where we have assumed 
that the densities pb — B 2 /(Sir) and p 7 are equal. In most circumstances (unless equipartition is very badly violated) 
Inverse Compton losses can be neglected. 

Charged pions (and kaons) traveling in a radiation field can also lose energy for photo-hadronic interactions. For 
completeness we have calculated the interaction time of pions in the GRB radiation field. The result is also shown in 
fig- El F° r this calculation the ttj hadronic cross section has been estimated as the sum of contributions for production 
of the relevant resonances [28[ : p(770), al(1260), 61(1235) and a2(1320) and a non resonant background. At high 
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energy (s > 30 GeV 2 ) the cross section is described by a formula obtained using Pomeron universality and Regge poles 
factorization from the fits to the 7-proton and 7r-proton cross sections reported in [281 ] . The lower energy non-resonant 
background has been normalized to obtain a smooth energy dependence. 

At very high energy pion-photon interactions are more important that Inverse Compton losses, however if the 
energy density in photons and magnetic field are comparable, the synchrotron losses dominate. 



A. Decay Energy Distribution 

To describe the effect of energy losses, it is useful to consider the function dpd cc /dE(E, Ei) that gives the probability 
density for an unstable particle of mass m and lifetime r, created with initial energy Ei to decay with energy E. 
Assuming that the energy loss is a continuous process, and can be described by the equation: 



dE 
~dt 



(E) = f(E) , 



(with f(E) is the average loss per unit time), the distribution dpdec/dE(E, Ei) can be calculated as: 



dE 



■{E-Ei) 



1 



exp 



dE' 



r J E f{E')E> 



r f(E)E 

For an energy loss of form f(E) — aE n (with n > 0) the decay probability becomes [2£ 

dpdec , 



dE 



■ (E; Ei 



E n+1 



with: 



exp 



m 
t a 



e 
n 



1 1 

E" ~ E? 



E/f(E) = l/(aE 



The quantity e has the physical meaning of the energy for which the loss-time: t\ c 
to the decay time: idecay = t E/m. 

For Ei -C e, the energy loss is negligible and the decay energy distribution becomes a simple delta function 



(37) 
(38) 

(39) 

(40) 
^ is equal 



dE 



■(E;Ei 



6[E-Ei], 



while for Ei 3> e the decay energy distribution takes a universal form 



dpd 



dE 



■(E; Ei 



E n+1 



exp 



1 e™ 
nE" 



(41) 



(42) 



An illustration of the decay energy distribution for the case f(E) = a E 2 , relevant for synchrotron emission is shown 
in fig. [7] The important qualitative feature in fig. [7] is that all particles with sufficiently high initial energy (Ei > 2 e) 
have nearly identical final energy distributions, with a well defined maximum at E ~ 0.6 e. In many cases this can 
result in "pile-up" effects in the final state particle spectra. 

VII. HADRONIC INTERACTIONS 

The target for the primary particles interactions can be either normal matter in gaseous form, or a radiation field. 
In the following we will consider both cases separately. 



A. Gas Target 

For a target material made of normal matter at rest the interaction rate of a primary particle of energy E p is: 

K pp (E p ) = n gas a pp (s). (43) 
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where n gas is the number density of the target material. Equation (|43[) is valid for a situation where the primary 
particles are protons and the target material is also dominated by protons. This is likely to be a good approximation 
in most circumstances. It is straightforward to consider the more general case. 

Since the target is at rest, the interaction rate is proportional to the cross section at a well defined cm. energy 
s = nip + 2m p E p . Since hadronic cross sections grow only logarithmically with cm. energy (or as power law oc s e 
with a small exponent e), the interaction rate of the primary particles changes only very slowly with energy. 

It is well known that the multiplicity and energy distributions of the particles produced in hadronic interactions 
cannot be calculated from first principles, however it is believed that particle production, to a a reasonably good 
approximation, satisfies Feynman scaling, defined by the condition: 

E d ^{jp\,Vs)=F j {x F ) (44) 
dp,, 

where pj, is the longitudinal momentum in the cm. frame, dtjj is the inclusive differential cross section for the 
production of particle type j, and xf = 2p*^/y/s is the Feynman variable. 

For large cm. energy (s 3> mj) the target rest frame energy of secondary particles in the forward hemisphere 
(xp > 0) is well approximated by the expression E ~ E p xf (where E (E p ) is the energy of the final state (projectile) 
particle in this frame). It follows that the validity of Feynman-scaling in the fragmentation region implies also the 
approximate validity of scaling of the inclusive cross sections in the target rest frame: 

^Ep { E; E p ) ~ ±- F pp ^ a (£) (for E » m) (45) 

The approximate validity of the scaling law (j45]l together with the slow variation of the hadronic cross sections 
with cm. energy have the important consequence (based on essentially the same argument previously discussed in 
connection with decay (|25[) ) that a power law spectrum of primary particles generates spectra of secondaries also 
having a power law form. For example for final state particles of type 7r: 

f 1 a-l 

Q-TriEiT^ C p Tlgas &pp E^ J dx X E pp — *7r{,^) — C p Tlg as G pp Z pp — ^(d) E n . (^6) 

Jo 

In general, one expects some deviations from a perfect power law. First of all, this can reflect a shape of the 
primary spectrum different from a power law. In addition the energy dependence of the cross sections, and the 
existence of violations of Feynman scaling (as measured in hadron colliders) introduce some distortions. Another 
situation that can result in large deviations from a power law emission arises when particles of different rigidities 
have different confinement volumes having different target densities (while our derivation of (|46| implicitely assumed 
a homogeneous source volume). 



B. Photoproduction 

In several of the proposed neutrino sources, the target of the primary particles is a radiation field. In this case 
the interaction rate and the energy distribution of the particles produced in an interaction depend not only on the 
density, but also on the energy and angular distribution of the target photons. The interaction probability per unit 
time of a proton of energy E p traveling in the radiation field described by n 7 (e, f2 7 ) can be calculated as: 

1 r r+ l dcosO 

K pl (E p ) = - — — — = de — -51 (1 - cos6> P7 ) n 7 (e, cos P7 ) a pl (e r ) (47) 

A pl yhi p ) J J _ 1 I 

where 8 pl is the angle between the photon and the proton momenta in the interaction, a pl is the photoproduction 
cross section, and e r is the photon energy in the proton rest frame: 

E £ 

e r = (l-cos6Ly) (48) 
m p 

The quantity e r is in one to one correspondence with the cm. energy of the reaction. It is convenient to change the 
integration variable from cos 9 P1 to e r . Restricting ourselves to a situation (and a frame) where the photon distribution 
is isotropic one can then rewrite (|4"T1) as: 

1 m 2 f°° f°° n (e) 

K pi( E p) = o~d / *r e r a pl {e r ) / de (49) 

A ^p Je th J(m p e r )/(2E P ) £ 
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where e t h is the threshold photon energy for pion production in the proton rest frame: 



ml 



eth = m, + ^ (50) 
Equation (|49|) can be recast in the form 

K P1 (E P ) = / de r fpy{e r ; E p ) . (51) 



This expression shows explicitely the fact that the pj interactions of a proton of energy E p do not correspond to a 
single value of the cm. energy but have a distribution, that in general is a function of E p , and is determined by the 
energy (and angular) distribution of the target photons. The probability distribution for e r is: 

A phenomenologically important case for the target radiation field is the form: n 7 (e) = CyE-P, that is an (isotropic) 
power law spectrum of slope (3. For this form the last integration in (I49p can be performed analytically with the result: 



9/9 / E \ l3 ~ 1 C°° 

K pl {E p ) = C 7 — - [-£-) / de r e-P a P7 (e r ) = K (J3) E^ 1 (53) 
p -t- 1 \m p / j eth 

One can see that the interaction rate has the energy dependence £^ _1 . The growth with energy (for /3 > 1) of the 
interaction rate can be understood observing that a proton of energy E p can interact inelastically only with photons 
above a minimum energy m p eth/ (2£7 P ). This minimum target energy decreases proportionally to E~ x , and therefore 
protons of higher energy can interact with a softer and more abundant photon population. A result essentially 
equivalent to equation (|53j) was originally shown by Waxman and Bahcall in [25| . An important feature of equation 
(|53"|) is that the probability distribution for e r takes a form that is independent from E p : 



p(e r ) =e r a pl {e r ) 



de r e r 13 cr p7 (e r ) 



i -l 



(54) 



This fact has interesting consequences discussed below. It should be noted however that equations (|53|) and (|54|) are 
calculated assuming a power law that extends to very large energy without any cutoff. This is not only unrealistic, 
but in some important cases also untenable. In fact the integral over e r in (|53|) diverges at its upper limit (since 
<7 P7 is slowly increasing with energy) for /3 < 1, and the expression for K pl becomes meaningless. The divergence 
corresponds to the divergence of the photon number density in the absence of an upper limit cutoff (note that in fact 
the energy density of the photon population diverges at the upper end already for j3 < 2). The introduction of a high 
energy cutoff for the energy distribution of the target photons is therefore mandatory. This will be discussed in the 
next section. 

The properties of particle production in jry interactions is clearly intimately related to the distribution of cm. energy 
(or equivalently e r ) of the interactions. It is necessary to consider the energy distributions of secondary particles in 
the "source frame" (where the primary proton has energy E p ). In this frame the energy of a secondary particle can 
obviously be expressed as function of quantities in the cm. frame with an appropriate Lorentz boost: 

E = 1 {E*+vp* z ) (55) 

where v and 7 are the velocity and gamma-factor of the cm. of the reaction in the source frame, E* is the cm. 
energy of the secondary particle and p* z the momentum component parallel to (3. The Lorentz 7 that connects the 
source and cm. frames is: 

E p + e E p 

7 = ^T"7i- (56) 

In the second approximated equality we have neglected the photon energy e with respect to E p , this is expected to 
be an excellent approximation. Similarly one can safely make the approximation v ~ 1. With these approximations 
we can rewrite the Lorentz boost from the cm. to the source frame as: 

E = E p ±J± =E P £. (57) 

V s 
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This equations indicates that to a good approximation all secondary particles of source frame energy E are created in 
the cm. frame of the interaction with (to a very good approximation) the same value of the quantity £ = (E* +p*)/ \fs. 
The energy spectrum of (for example) pions created in the interaction of a proton of energy E p can then be written 
as: 



^P^(E V ; E p )^± rde r p(e r] E p 



dn 



P7 — >7T 



(58) 



i=E w /E p 



where dn p7 ^ 7r /<i^(^, e r ) is the £ distributions of the secondary particles of type n in fry interactions with cm. energy 
that correspond to e r ; this distribution is convoluted for a fixed value £ = E v / ' E p over all possible value of the e r with 
the appropriate distribution. If the energy distribution of the target photon field has a power law form, the function 
p(e r ; E p ) is independent from E p , and equation (|58p becomes the expression of a scaling law of form: 



dn 



dE„ 



(E n ; E p ) 



E n Fp1 ^ { E ; 



E^ 



(59) 



The scaling function Fp-y—,^ is not truly universal, but depends on the slope of the target photon spectrum. Integrating 
over all primary particle energies one can obtain the production rate of pions as: 



Qtt(Ett) 



dE p Np(Ep) K pl {Ep) 



dn r 



dE-w 



-(En-; E p 



dEp (C p E- a ) (K (p) Eg- 1 ) 



E„ 



F, 



C p K {(i)E-^-^ / d£ C' P Fp^it 0) 
Jo 

CpKo(P) Z pl ^[a-P +!,/?] E-^-P+V . 



(60) 



This equation shows that under the assumptions made, that is: (i) a power law spectrum of protons, and (ii) an 
isotropic power law spectrum of target photons, the energy distribution of produced secondaries is again a power law 
with slope a — j3 + 1. Note that this result is based on purely kinematical considerations, while the result (|46|) about 
the interactions on an ordinary matter target was based on a dynamical assumption about the approximate validity 
of Feynman scaling in hadronic interactions, and on the weak energy dependence of hadronic cross sections. 

A result roughly equivalent to ([50|) was obtained by Waxman and Bahcall in [25|, who concluded that in the case 
[3=1 the neutrino emission has a power law spectrum with the same slope of the parent protons. However, the 
validity of equation (|60|) . as the validity of equation (|53[) . relies on the untenable assumption that the power law 
spectrum of target photons has no high energy cutoff. The important effects of the existence of an high energy cutoff 
for equations (|5"3")1 and (pO)) will be discussed in section [Villi 



C. Nuclear photodisintegration 

The beta decay of neutrons produces T> e . In many circumstances the contribution of this anti-neutrino source is of 
negligible importance. This is because the average fraction of the neutron energy E n carried away by the V e (in a frame 
where the neutron is ultrarelativistic) is (Ex> e )/E n ~ 5.1 x 10 -4 , with end-point £;™ ax / E n ~ 2(m n — m p — m e )/ra n ~ 
1.66 x 10 -3 . In most cases the flux of the softer rt-decay neutrinos is neglible with respect to the contribution of 
neutrinos from tt/K of decay that carry a much larger fraction (of order ~ 0.25) of their parent energy. There are 
two circumstances where the n contribution can become important. The first one is when the neutrino spectrum has 
a low energy cutoff. This happens naturally when the target of the primary particles is a radiation field and there is 
an interaction energy threshold. In this case the v e produced in neutron decay can become the dominant component 
of the neutrino flux at low energy, because most of the neutrinos from pion and kaon decay have higher energy. An 
example of this situation will be shown in the following section. A second, more interesting case is when the neutrons 
are produced in the photodisintegration of high energy nuclei [HJ . In this case it is in principle possible to have the 
emission of a pure V e flux. In fact the threshold for photodisintegration of a nucleus of mass number A expressed in 
terms of energy per nuclcon (E = E tot /A) is of order: 

{EoV A ^ r^bind (61) 
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where £ 7 is the energy of the target photons and ebind — 8 MeV is the binding energy of a nucleon in the nucleus. 
The threshold for pion photoproduction is: 

Since the binding energy ebind is approximately fifteen times smaller than a pion mass, it is in principle possible to have 
circumstances where the primary particles are below the threshold for pion production, but above the threshold for 
photodisintegration. This would clearly result in a pure V e flux. It should be noted that the photodisintegration and 
pion photoproduction thresholds differ by only one order of magnitude, and therefore a pure V e flux can only extend 
for a small interval of energy. In general one expects important contributions from neutrinos from meson decays. 
In fact, even a very small number of pion production interactions created by the high energy tails of the photon 
and/or primary particles spectra can result in an important "contamination" of v e and muon neutrinos. Most of these 
neutrinos have an energy two orders of magnitude larger that the V e produced in n decay and have a correspondingly 
larger cross section. 



VIII. NEUTRINO PRODUCTION IN GAMMA RAY BURSTS 



As an illustration of the spectra and flavor composition of neutrinos from astrophysical sources, in this section we 
will consider the model of neutrino production in Gamma Ray Bursts (GRB) developed by Waxman and Bahcall 
(WB) [25]. Our goal is not to analyze critically the assumptions that underline the model, but to recalculate with a 
more detailed modeling of particle production the neutrino fluxes starting from the same general assumptions. 

In the (WB) model [25| the GRB prompt emission is generated by the synchrotron radiation of high energy electrons 
accelerated by internal shocks in a relativistic expanding wind (for more discussion see [33|,[34[; alternative explanations 
of the GRB mechanism exist in the literature, see for example [35|). 

The internal shocks also accelerate protons, that can interact with the radiation field of the wind and produce 
secondary particles that generate neutrinos. The accelerated protons have a power law spectrum of slope a ~ 2, while 
the target radiation field is isotropic in the wind frame, and approximated as a broken power law of form: 



n 7 (e) 



C-y £ 

C~Y £\) £ 



for e < £b, 
for e > £b 



(63) 



where £b is a break energy. The photon energy distribution in the wind frame is related to the obervable spectrum of 
the GRB prompt emission by a Lorentz transformation with parameter T ~ 300. Since the observed break energies 
are distributed around an average value (£b bs ) ~ 300 KeV [36j, the typical value of £b in the wind frame is of order 
1 KeV. 

The energy spectra of individual GRB are indeed well fitted by two power laws smoothly joined at a break energy, 
however the fitted values of the exponents of the s pec trum have rather broad distributions centered at ~ 1 below the 
break energy, and (3 ~ 2.2 above the break energy [361 ]. In our calculation we have therefore considered a generalization 
of ([(521) that leaves the two exponents as free parameters. We have also introduced a high energy cutoff £ max , because 
in its absence the energy density of the radiation field would diverge for a high energy slope < 2. Our description of 
the energy spectrum of the target radiation field then becomes: 



(C 7 /e b ) . 
(C 7 /e b ) (e/et)-* 




for e < £ b , 
for £ b < £ < £ n 
for £ > £ max . 



(64) 



The calculation of the interaction rate of a proton of energy E p in such a radiation field is straightforward. The target 
photon field (|64[) can be seen as the sum of two power law spectra with sharp low and high-energy cutoffs £ m i n and 



We have already obtained the expression for the proton interaction rate for a power law that extends to all energies. 
The expression ([53| must however be modified in the presence of low and high energy cutoffs. In the presence of 
a high energy cutoff £ m ax the interaction rate vanishes below a threshold energy Eth- Above the threshold, in the 
interval: 



Eth — 



2f ~ p ~ If ■ 

^ c max ^ c min 



the interaction rate is 
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(65) 



(66) 
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coinciding with (153)) in the limit e max — ► oo. The effect of a lower energy cutoff e m i n is to stop the E@ 1 growth of 
the interaction rate. For E p > m p eth/(2£ m i n ) equation (|66p must be replaced by: 



Cry Tn^ 



K„{E p ) - 2 



(1 + 13) 
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(67) 



For a qualitative understanding, it can be useful to consider the ^7 cross section as approximately constant above 
threshold. The integrations over e r in (f6"6")) and (|6"Tj) are then trivial. For j3 ^ 1 one finds: 
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(68) 



where x = E p / E t h is the proton energy expressed in units of the threshold energy, and r 
finds: 

for x < 1 

[K PJ (E p )f =1 ~ <| C 7 CTp7 [log[x] -|(1- £)] for 1 < x < r 



for x > r 

£max/£min- For /3 = 1, one 



C 1 (T pl [k>g[ 



(69) 



for 2! > 



These expressions show that the interaction rate for a power law spectrum of target photons in a restricted range 
do maintain the approximate behaviour K pl oc E^ 1 but only in a limited energy region (in the limit /3 — ► 1 the 
interaction rate grows logarithmically: _Kp 7 oc log[_E p e max ]). The interaction probability vanishes below threshold 
and goes asymptotically to a constant for very large E p . It can also be useful to consider the high energy limit 
(E p /Eth — > oo) of expressions (|6"8|) and (f6"9"| : 
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(70) 



[tf P7 (£ p )] 



t P7 C 7 log 



(71) 



These expressions have the form K pl — a pl N 1 with iV 7 



the integrated number density of the target photons, that 
can be immediately recognized as the correct high energy limit. 

The interaction probability for the radiation field (|64[) can be obtained combining two expressions corresponding 
to the parts of the photon spectrum below and above the break energy £b that plays in the two cases the role of the 
maximum or the minimum target photon energy. An illustration of the energy dependence of the proton interaction 
rate is shown in figure The curves in this and the following figures are made independent from the value of £b 
measuring all energies in units of: 
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The energy E* has the physical meaning of the threshold proton energy for inelastic interactions with photons having 
the break energy e\,. In figure [5] we have chosen for the target photons energy spectrum the original form (|6"3")) with 
exponents /?i = 1, (3% = 2 and no high energy cutoff. For low E p the interaction rate grows approximately linearly with 
energy according to the scaling law (|53f) because low energy protons interact inelastically only with the high energy 
part of the radiation field (with slope fa = 2). For larger energy (E p > E*) the growth becomes logarithmic following 
the qualitative behaviour of equation (|6l)|) . It should be noted that the assumption made in 25] of considering the 
interaction rate as approximately constant for E p > E* is not a good approximation. 

In fig. [5] we show the distribution f(e r ; E p ) of the cm. energy of the interaction for different values of the proton 
energy (E p /E* — 1, 10, 100). The distribution is proportional to theory cross section, and the shape of the curves in the 
figure reflects the energy dependence of a pl that has prominent resonances. The most important one (at e r ~ 0.3 GeV) 
corresponds to the production of the A resonance. The important feature in the figure is that with increasing energy, 
higher and higher cm. energies become possible. This has significant phenomenological consequences. For example, 
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for e r sufficiently small, only production of a single pion, in the two channels pir° and nir + , is kinematically allowed, 
and therefore there is no ir~ production, and consequently no V e production. For larger cm. energy multiple pion 
production becomes possible, V e can be produced, and the ratio v e /v e decreases. Above the KA energy threshold, 
also the production of kaons becomes possible introducing an additional neutrino source. 

In the following we will show some examples of the neutrino fluxes that are obtained varying the parameters of the 
model. Our calculations have been performed by Montecarlo methods using a detailed model for the ^7 cross section 
that includes the production of all relevant resonances and a non-resonant component [381 ]. In order to perform a 
calculation one needs to specify: (i) the primary proton spectrum, (ii) the target photon spectrum, and finally the 
value of the magnetic field B (all other sources of energy losses for secondary particles are considered as negligible). 
The target photon energy distribution is taken with the form (I64[) that depends on 4 parameters: {C 7 , £b, e m ax, Pi, 
P-i\. Different values of £b correspond to different values for E* oc e^ 1 . The proton spectrum is taken with the form: 



that is a power law with a smooth high energy cutoff: defined by the three parameters: {C p , a, -Emax}- In all of the 
following we have kept fixed two of the parameters, the proton maximum energy: E max /E* = 10 4 , and the photon 
high energy cutoff e max /eb = 300. The parameters C p and C 7 determine the absolute value of the neutrino emission, 
but not its shape, and will be left unspecified in the following. 

The absolute value of the neutrino flux from an individual source can be obtained specifying the values of C p (that 
fixes the amount of energy in relativistic protons in the wind frame), C 7 that specifies the density of the photon 
target field, the exact Lorentz Boost and obviously also the source redshift (distance). The integration over the 
ensemble of all sources requires additional assumptions about the distributions of the relevant parameters used for 
their description and of their cosmological evolution. In this work we will not discuss this integration. 

The value B of the magnetic field is in a one to one correspondence with the muon synchrotron energy, defined in 
equation (|33j) . that can also be expressed as: 



As an intermediate step toward the calculation of the neutrino fluxes we show in fig. [TO] the yields of different 
mesons. The meson yields are calculated for a proton spectrum with exponent a = 2, and a radiation field with slopes 
(below and above the break energy) Pi = 1 and P% = 1. The solid lines in fig. [TO] show the energy distributions of the 
different mesons at the moment of their creation. In the absence of significant energy losses the curves also describe 
the energy distributions at decay. One can see that all possible mesons are produced, but with different abundances. 
The most abundant particles are charged pions, reflecting the fact that a large fraction of the pr/ interactions happens 
close to threshold where only single (neutral or positive) pion production is possible. The production of ir~ becomes 
possible only above the threshold for two pion production, and is therefore suppressed. The production of kaons 
is significantly smaller, because it is suppressed both dynamically (since it requires the creation of ss pairs) and 
kinematically (because of the larger strange hadron masses) . The threshold for the production of mesons containing 
a strange antiquark s (K + and K°) is lower, since it corresponds to final states containing a strange baryon (such as 
AK + ). The production of mesons containing a strange quark s (K~ and K ) has a higher threshold since the final 
state must contain a minimum of two kaons (as for example: pK + K~). Accordingly one finds that the productions 
of different kaon types is ordered as follows: K + > Kl > K~ . The shapes of energy spectra of the mesons are never 
well approximated by a simple power law, but are always "curving" in a log-log representation. The shapes of these 
spectra reflect the energy dependence of the proton interaction rate, and the opening up of the different kinematical 
channels. At the highest energy the rapid drop in the meson yields is connected to the cutoff in the primary proton 
energy at E max /E* ~ 10 4 . Note that the ratio tt~ /ir + and K/ir are not constant but increase with energy reflecting 
the increasing importance of 7r~ and kaon production with growing cm. energy. 

The dashed lines in fig. [TO] describe the decay energy distributions of the produced mesons assuming the presence 
of a magnetic field that corresponds to e M = 3 (or B = 2.8 x 10 3 £^ cV Gauss). Pions created with energy above 
their critical energy for synchrotron losses (in units of E*\ e n ~ 18.4e M ~ 55) lose most of their energy before decay, 
and therefore the number of pions decaying above this critical energy is strongly suppressed (oc (E v / e^) -2 ) ■ One 
can also notice that synchrotron losses result in an enhancement of the number of pions decaying with energy just 
below €„. This is a pile-up effect due to the fact that all high energy pions (with E v ^> e n ) decay with similar energy 
distributions. The same effects are present also for charged kaons at higher energy, because the synchrotron energy 
for charged kaons (in unit of E*) is ck — 628e M ~ 346^ ~ 1900. 

Secondary neutrons are also produced in pj interactions. Assuming that they can freely exit from the source, the 
V e spectra from their decay have been calculated and included in the following figures. 
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Some examples of the resulting neutrino fluxes obtained summing over all possible parent particles and all v types 
are shown in fig. [11] and fig. [T3J The lines in fig. QT] are calculated assuming that energy losses for secondary particles 
are negligible. The different curves correspond to different assumptions about the shape of the energy spectra for the 
primary protons and target photons. The thick solid curve is calculated with the choice of slopes (for the proton flux 
and for the radiation field below and above the break energy): {a, (3i, /3a} = {2, 1, 2}. The other three curves show the 
effect of changing the slopes one by one. The largest effect is related to the slope of the primary protons, increasing 
a from 2 to 2.4 results in an important softening of the neutrino spectrum. The modification of the shape of the 
energy distribution of the target photons also distorts the neutrino spectrum. Changing the slope below (above) the 
break energy, modifies the high (low) energy part of the neutrino spectrum. This is easily understood, since high 
(low) energy neutrinos are produced by higher (lower) energy protons that mostly interact with lower (higher) energy 
photons. Note that (reflecting the spectra of the parent pions and kaons) the neutrino energy distribution changes 
gradually its slope and cannot be well fitted by a power law. 

Figure [T2l illustrates the effect of including synchrotron losses on the neutrino fluxes. The different lines in fig. [T2l 
show the neutrino flux (summed over all v types) calculated with the slopes {a, 01,02} = {2,1,2}, and three different 
assumptions for the magnetic field that correspond to muon synchrotron energy e M = oo, 30 and 3 (the first case 
corresponds to negligible losses). Increasing the magnetic field (reducing e^) suppresses the neutrino flux at high 
energy. In general, in the presence of a strong magnetic field in the source, one can identify four interesting neutrino 
energy ranges, that are related to the ordering of the synchrotron critical energies for different particles: e M < e ff < ex- 
At sufficiently low energies the synchrotron losses are completely negligible; at higher energy the muons energy losses 
have to be taken into account; at still higher energy also the losses of charged pions must be considered; at the 
highest energy the emission of synchrotron radiation is important also for charged kaons. In general, in the presence 
of important synchrotron emission, one can have some pile-up effects, due to the fact that all high energy particles 
decay just below their critical synchrotron energy. 

The inclusive neutrino spectra shown in the previous figures contain the contribution of the v e from n decay. The 
structure of this contribution is illustrated in fig. 1131 The anti-neutrinos from neutron decay have significantly lower 
energies that the neutrinos from meson decay. They are the main component of the neutrino flux at the lowest 
energies, where however they are difficult to observe because the flux is suppressed and the cross section is small. 

The flavor ratios R efl = (u e +F e )/ (y^ +^ M ) that correspond to the neutrino spectra of the previous figures are shown 
in figures [T4l and [T51 The flavor ratio is energy and model dependent. The i? e/J ratios shown in figure [T^l correspond 
the energy spectra of fig. [11] and are calculated assuming that all energy losses are negligible. At the lowest energy 
the ratio R ell increases rapidly due to the contribution of V e from n decay. At higher energy the ratio i? eM stays close 
to the value of 1/2, however the expanded scale allows to see that the "naive" result is not exact. The value of i? e/i , 
for the reasons illustrated in section [V] are correlated with the shape of the energy spectrum. The gradual softening 
of the neutrino energy distribution with growing E v is reflected in a slow increase of i? eM . Comparing the different 
curves, one can also see that the curve calculated for a — 2.4 that gives the softest neutrino spectrum correspond (at 
sufficiently high energy) to the highest i? eM . 

The effects of the existence of significant synchrotron losses are very important for the flavor ratio. The main effect 
is the effective loss of the neutrinos from high energy muon decay. This results in a "pile-up" effects that produced 
an increase in the ratio for neutrino energies below the synchrotron muon energy e^. Above this energy the ratio 
falls steeply reflecting the effective absence of the main source (muon decay) of electron neutrinos. The flavor ratio 
however does not vanish because kaons have decay modes into electron neutrinos, and reaches a level of a few percent. 
At very large energy the role of neutral kaons (that do not suffer synchrotron losses) become enhanced and the flavor 
ratio grows because of the importance of the K ° 3 decay mode. 

Figure \W\ and figure Q~7] show the neutrino/ anti-neutrino ratios v e jv e and v^jv^ calculated for the same three 
situations present in fig. W2\ and [T51 (that is slopes {a, 0i, fa} = {2,1,2} and three values of the magnetic field that 
corresponds to e M = oo, 30 and 3). 

The vjv ratios are difficult to measure because the detectors have no capability to measure the charge of the final 
state charged lepton in charged current neutrino interactions. The most attractive idea [53] to estimate the v e /v e ratio 
is the comparison of the event rate at and near the "Glashow resonance" [3!| at E v ~ M^ v /(2m e ). The resonance is 
present only for V e in reactions such as V e + e~ — ► hadrons that can proceed via the formation of a W~ boson in the 
s channel, while non-resonant events proceed via the normal charged current interaction on nucleons. 

The study of v e /v e ratio is particularly interesting because it has been proposed [37] as a method to distinguish pp 
and p7 interactions as the source of the neutrinos. The idea is that the v e /v e ratio at the source is close to unity for 
a pp source, when the interactions produce approximately equal numbers of positive and negative pions (that decay 
into v e and z7 e ), and is very large (or actually diverges) for pry interactions, when most interactions happen close to 
threshold, and the cross section is dominated by single pion production (via the A resonance), and therefore only 7r + 
production is present. This argument is qualitatively right, but quantitatively incorrect, in fact as discussed above, in 
general increasing E p a broader range of cm. energies becomes possible, and this results in a growing contribution of 
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negative pions to neutrino production. In certain circumstances also the contribution of kaons can become important. 
This is illustrated in fig. rTBjthat shows the v e /v e ratio. At the lowest energies the flux is dominated by the v e from n 
decay, and the v e /v e ratio is therefore small. When the contribution of neutrinos from pion decay becomes dominant 
the ratio v e /v e is significantly larger than one, because positive pions are more abundantly produced than negative 
ones. In fact below the threshold for two pion production the ratio n + /n~ diverges. With growing energy the ratio 
7r + /n~ decreases because of the growing importance of multiple pion production. This is reflected in a decrease of the 
v e /v e ratio that at large energies becomes ~ 2. If a large magnetic field is present, at large energy the synchrotron 
energy losses for muons become important and neutrinos from muon decay are suppressed. The decay of neutral 
kaons (Kl) becomes the dominant source of v e and v e , and the ratio v e jv e becomes unity. 

The ratio v^jv^ (shown in fig. [T7|) is approximately unity at low energy, reflecting the fact that each charged pion, 
after chain decay, contributes a and a V^, The existence of important synchrotron losses can modify this quite 
robust prediction. At high energy, when muons energy losses become significant, only the neutrinos from direct pion 
decay are important, and the ratio v^jv^ is approximately equal to the tt + /tt~ . At larger energy, charged kaons 
become the dominant and source and the v^jv^ ratio increases reflecting the higher K + /K~ ratio. 

In fig. [15] we show the {v e /v^) ratio for an Earth observer. To propagate neutrino we have used standard oscillations 
with mixing parameters #23 = 45°, 612 = 34° and 6*13 = 0. In the figure we have also shown the same observable ratio 
for a source with the "naive" composition [v e . v^, v T \ = [1, 2, 0]. For the "best fit" values of the mixing parameters the 
observable (i/ e /^ M ) has the well known value of unity. Changing the value of the 6*23 angle by ±5° the ratio changes 
by —13% (+7%). Similarly for 6*13 = 5° the ratio take values in the interval (0.95,1.04) depending on the value of 
the phase 6. This figure illustrates the interplay between astrophysical and particle physics uncertainties (also shown 
before in equations ([3]), (T5|) and ([5])). In order to perform a measurement of 623 one clearly needs an accurate control 
of the initial flux composition. 

The description of the emission of neutrinos from GRB that we have used in this section following reference [25j 
is very simplified, and could be improved in many aspects. The results of our calculation show that one obtains non 
trivial and intriguing energy dependences of the flavor ratios, that could reveal important informations about the 
source. This qualitative observation is of a much more general validity and one expects that in most cases, and also 
using more realistic descriptions of the sources, one should find non trivial structures in the energy spectra and flavor 
composition of the astrophysical neutrino signals. 

IX. EXPERIMENTAL FLAVOR IDENTIFICATION 

In this work we will not consider in detail the crucial problem of the experimental determination of the flavor ratios. 
Several experimental approaches for the detection are at present being developed. The concept that is in the most 
advanced stage is a large volume (cubic Kilometer) ice or water Cherenkov detector, with photon detectors distributed 
inside the volume, as originally proposed b y th e Dumand group [4fj| and now in different stages of development at 
the South Pole [U, \^ , in the Baikal lake [43l | and in the mediterranean sea [44l |45| . |4(| . Several other alternative 
techiques are also being developed that includes radio [H,|4?|, acoustic and air shower [5(J [EU detection methods. 

A very imporant problem is of course the amount of data that is required in order to reduce the statistical errors 
to the level needed for a meaningful measurement. Today we only have upper limits for the astrophysical neutrino 
fluxes, optimistically one can expect that these fluxes will be soon discovered just below the current upper limits, 
however even in this case the event rates will remain small, and one will need very long data taking periods, or much 
larger mass detectors to collect sufficient data [H, HH • 

In this section we want to stress the point that in all the experimental methods that are envisaged, the measurement 
of the flavor ratio is connected with a determination of (or a theoretical assumption on) the shape of the energy 
distribution of the neutrino fluxes. This happens because the detection methods of different neutrino flavors have 
efficiencies with different dependences on E v . If the statistical errors are sufficiently small (with very large mass 
detectors), and having achieved perfect control of the detector acceptances and efficiencies, the uncertainties in the 
determination of the neutrino spectral shapes will remain as the dominant source of error in the measurement of the 
flavor ratios. 

To illustrate the interplay between the spectral shapes and the determination of the flavor ratio one can consider 
the case of the Km 3 type Cherenkov telescopes. In these detectors one expects to find different classes of neutrino 
induced events [l3| . 

(i) In "Track Events" a single up-going muon is seen entering the detector. The energy of the muons can be 
estimated from the amount of Cherenkov light produced by the photons and pairs radiated from the muon 
track. 

(ii) In "Shower Events" one observes the release of a large amount of Cherenkov light generated by the shower 
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induced by a neutrino interaction in the detector volume. The amount of light can be translated in the estimate 
of a visible energy release i? v i s . 

(iii) In "Double Bang Events" [1, HiJ one detects two distinct energy releases inside the detector volume, with space 
and time separations that are consistent with the propagation and decay of an energetic tau lepton. 

Events in class (i) are associated with the charged current (CC) interactions of muon neutrinos below the detector. 
Events in class (ii) can be produced by several types of neutrino interactions. Electron (anti)neutrino CC interactions 
are the main source of shower events, there are also contributions from the neutral current interactions of all neutrino 
types, and in general also from the CC interactions of and v T . Finally events in class (iii) are generated by v T 
CC interactions. Clearly, the ratio of the frequencies of these three classes of events does give information about the 
flavor composition of the neutrino fluxes. In fact the ratio of the event rates for "Shower" and "Tracks" events has 
been proposed as the best method for a measurement of the v^jv^ flavor ratio. 

The frequency of "Track" events can be studied as a function of the visible muon energy E^. The differential event 
rate can be calculated as: 

dE N dT S dt {EM = A Sr k8 (^^e) E / dE « ME.M E v ) (75) 

where A^,^ a (E IM) f2®) is the detector effective area for the detection of muons with energy E^ from the direction 
Q®, dY^^ fl /dE fl (E fl - E v ) is the "muon yield", that is the probability that a (z7 M ) of energy E v produces a fi T of 
energy E^ at the detector. In this equation we have assumed that the muon and the neutrino are collinear. Even 
for an individual source, the energy spectrum will depend on the zenith angle (9®, because of the effects of neutrino 
absorption in the Earth. This procedure hower implies a good understanding of the neutrino energy spectrum. 
The rate of "Shower" events with a visible energy release E v - ls can be estimated as the convolution 
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In this equation Mdet(-E'vis) is the effective detector mass; the first line is the contribution from the charged current 
interactions of electron (anti)-neutrinos, in this case the visible energy, in first approximation (neglecting differences 
in the light to energy conversion for cm. and hadronic showers) equals the initial v energy; the second line is the 
contribution from neutral current interactions, in this case the visible energy is only a fraction of the initial energy 
because the neutrino in the final state carries away an invisible energy E v (1 — y). The contributions of and v r CC 
interactions to the shower event rate in (|75|) has been left implicit. A fraction of these interactions can probably be 
identified and eliminated selecting events with a muon in the final state, or a double bang structure, but the rest will 
contribute to the "shower" class of events. 

From ([75]) and (|76p it is clear that the energy response for the two classes of events is different, with the "Track" 
events selecting larger E v . Therefore the ratio for the integral rates for "Shower" and "Tracks" events is strongly 
dependent on the spectral shape of the fluxes. Use of the differential distributions in E^ and E Y - m allows in principle 
to unfold the shape of the neutrino energy spectra, however an important difficulty is that a narrow interval in the 
muon energy E^ (that can however be measured only with modest resolution) corresponds to a very broad interval in 
E v . Also E'vis in shower events is not in a one to one correspondence with E v . The unfolding procedure has therefore 
some significant limitations. 

As an illustration of these problems, one can look back at the original discovery of neutrino oscillations [f| 0, [5f| 
with atmospheric neutrinos. The evidence collected by the SK [55[ experiment for subGeV fully contained events 
(-Evas < 1-33 GeV) implied comparing the ratio of the observed rates of e-like and /j— like events with a MC prediction 
(A*/e)MC — 1-50. This prediction is significantly smaller than the no-oscillation ratio (y^ + V^)/(v e + V e ) ~ 2 for the 
combination of several reasons, related to the difference in detector acceptances for e-like and ^-like events. The 
experimental cuts select for the e-like events a broader energy range, both at the low and high energy ends (at low 
energy because electrons have a larger yield of Cherenkov photons, and at high energy because electrons have a shorter 
range and are more easily contained). For this reason the MC prediction does depend on the shape of the neutrino 
fluxes. For these sub-GeV events, the intervals of E v for the two classes of e-like and //-like events are very similar, 
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and therefore the associated systematic uncertainty in the (/i/e) prediction is relatively small, but it remains as a 
relevant source of error. 

The comparison of the rates of classes of events that correspond to very different atmospheric neutrino energy, such 
as samples of fully contained e-like events and upgoing muons (that are produced by charged current interaction 
with a median energy of order 100 GeV) also allows in principle to verify the existence of <-> v T transitions. However 
in this case the systematic uncertainty in the prediction of the shape of the energy distribution is much larger, and 
in practice this method gives only marginal evidence for the existence of a suppression of the upgoing muon rate. 

In conclusion, the extraction of the e/fi flavor ratio from the comparison of the rates of "Track" and "Shower" 
events in high energy neutrino telescopes requires excellent control of the detector performances, and a sufficiently 
precise knowledge of the shapes of the neutrino spectra. The determination of these shapes could be obtained from 
the data themselves (the "unfolding method" ) or with the help of reliable theoretical models for the source. In both 
cases the task appears as remarkably difficult. 

X. CONCLUSIONS 

The study of astrophysical neutrinos can be used for two distinct goals: the understanding of their sources and 
the investigations of the neutrino properties. These two goals are in a sense contradictory. The investigation of the 
neutrino fundamental properties requires a comparison of the observations with expectations that must necessarily be 
based on a sufficiently accurate understanding of the source. On the other hand the flavor composition of a neutrino 
signal encodes information about the source properties, but in order to extract this information one needs to know 
the properties of flavor transitions over galactic or cosmological distances with sufficient accuracy. 

Several recent works [ill, EH EH have speculated that it is possible to discovery new neutrino properties, or even 
measure precisely the parameters of the already established standard oscillations [56[ from the observation of future 
astrophysical neutrino sources assuming that the flavor ratios at the source are robustly predictable. The often 
repeated argument is that a "standard neutrino source" is dominated by charged pion decay, and therefore emits 
the different flavors (summing over v and v) with relative intensities: [v e ,i' fl ,i/ T ] ~ [1,2,0]. This statement however 
must be understood as a first approximation of only limited validity. First of all it requires that the energy losses of 
the muons created in charged pion decay are negligible. This crucial assumption is satisfied for several but not all 
of the proposed neutrino sources. Even for "thin environment sources" (where the density of ordinary matter and 
magnetic field is sufficiently low so that the energy losses for all secondary unstable particles can be safely neglected) 
the flavor ratios are (slowly varying) functions of the neutrino energy, with values that are determined by the shape of 
the neutrino spectra and by the nature of target material. This dependence arises from two effects: the first is that in 
general the contributions of other weakly decaying particles (kaons and in some cases also neutrons) is not negligible, 
the second is that, even if charged pions are the only significant neutrino source, the three neutrinos that are emitted 
in a 71"^ chain decay have different energy spectra. Folding these decay spectra with the energy distribution of the 
parent particles, one obtains slightly different shapes, and therefore energy dependent flavor ratios. In the special case 
where the parent pion spectrum is a power law, the neutrino energy distributions have also a power law form with 
the same slope, and the v flavor ratios are constants. As a (phenomenologically important) example, the neutrino 
fluxes with spectrum cx E~ 2 emerging from a "thin" , charged pion source, have relative intensity (at the source) 
[v e , i/a, v T \ ~ [1, 1.86, 0], a 7% deviation from the naive expectation. 

In general, the combined effects of different shapes for the v energy spectra, and different contributions of the relevant 
weakly decaying particles correspond to variations of order 10% of the Ve/v^ ratio (in some special circumstances 
even larger). The size of these corrections is however comparable to the effects on the observable flavor ratios due to 
the present uncertainties in the values of the standard oscillation parameters (for a quantitative illustration of this 
point see equations ([3]), Q and |5])). The very ambitious program [56| to attempt a precision measurement of #23 
or a determination of #13 with astrophysical neutrinos cannot therefore ignore these problems, and requires a very 
precise control of the structure and properties of the neutrino sources. 

Much larger effects on the neutrino flavor ratio arise from situations where the environment of the source contains 
large energy densities in magnetic field, or in ordinary matter. In the presence of sufficiently large magnetic fields, 
such as those that could exist in a GRB fireball, high energy muons can lose most of their energy to synchrotron 
radiation (29l . l3lj and their contribution to the neutrino fluxes becomes strongly suppressed. This results in a very 
small Ve/v^ ratio, with a value determined by the size of the kaon contribution (that generate electron neutrinos in 
direct decay modes). At higher energy also pions suffer important synchrotron losses, and kaon decay become the 
dominant v source for all flavors. In some extreme circumstances neutral kaon decay will remain as the only significant 
neutrino source with a v e jv^ ratio at the source of order unity. 

In some circumstances it is possible that interactions with ordinary matter play a major role. If the density is 
sufficiently large high energy pions reinteract before decaying. This supresses the neutrino fluxes and enhances the 
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contribution of kaon decay. Summing over charged and neutral kaons, this also results in a significantly reduced v e l v y. 
ratio. 

The prediction of the neutrino flavor ratios in these circumstances (strong magnetic field or large densities of 
ordinary matter) is going to be very difficult, because the assumption of a "single volume" homogeneous source is 
very likely a poor (and in any case not demonstrable) assumption. The neutrinos will likely emerge with different 
spectra and flavor ratios from different regions of the source, and for example fitting the data including a value of the 
field B, and/or a matter density n gas is likely to be only a poor approximation. 

The existence of neutrino sources dominated by neutron decay are also a theoretical possibility, in situations where 
photodisintegration of relativistic nuclei interacting on a radiation field is possible while pion photoproduction is 
kinematically forbidden. If such sources exist, it will however be very difficult to use these fluxes for studies of 
the neutrino flavor transitions because in general one expect them to have an important "contamination" of muon 
neutrinos in their initial flavor composition. In fact it is unnatural to have the pion photoproduction to be absolutely 
forbidden, because the energy thresholds for nuclear photodisintegration and pion production differ only by one order 
of magnitude. The size of the contribution of pion decay to the neutrino fluxes in these circumstances is difficult to 
control. 

Even considering all the caveats of our discussion, there are some circumstances where observations of astrophysical 
neutrinos could lead to the discovery that new physics is involved in the pro pag ation of neutrinos over astrophysical 
distances. In particular the investigation of the existence of neutrino decay [lq is a very attractive possibility. The 
observation of the absence (or very strong suppression) of a v e component in the measured fluxes cannot be explained 
with even extreme variations of the source astrophysical properties, and is on the contrary predicted (for any type 
of source) in a neutrino decay model with inverse mass hierarchy. For the normal mass hierarchy the neutrino decay 
model predicts a large {y e / v ii)ohs ratio between 3 and 6. If the neutrinos have standard properties, ratios (ve/^^obs 
much larger than unity are only possible (and limited to < 4) for sources dominated by v e emission. Therefore from 
the observation of a large (i / e /^ A ,) bs ratio it could be possible to infer the existence of v decay in the normal mass 
hierarchy case. This could require some understanding of the properties of the v source. 

In the general case the observation of flavor ratios very different from the "standard source" values (relative inten- 
sities [v ei Ufj,, v T ] bs ~ [1, 1, 1]) could be interpreted as evidence for some form of new physics in neutrino propagation, 
but also, and probably more economically, as an indication of the existence of some unexpected properties in the v 
source. Resolving these ambiguities is in principle possible, but certainly not easy, and requires a careful interpretation 
of the multi-messenger (photons, neutrinos and possibly cosmic rays) and multi-wavelength observations. 



[1] T. K. Gaisser, F. Halzen and T. Stanev, Phys. Rept. 258, 173 (1995) [Erratum-ibid. 271, 355 (1996)]; 
[2] J. G. Learned and K. Mannheim, Ann. Rev. Nucl. Part. Sci. 50, 679 (2000); F. Halzen and D. Hooper, Rept. Prog. Phys. 
65, 1025 (2002); 

[3] P. Lipari, Nucl. Instrum. Meth. A 567, 405 (2006) [astro-ph/0605535| . 
[4] EGRET Collaboration, Astrophys. J. Suppl. 12 3 79 (1999). 

[5] HESS Collaboratio n, Science 307, 1938 (2005) [astro-ph/0504380| . HESS Collaboration, Astrophys. J. 636, 777 (2006) 
[astro-ph/0510397] . 

[6] K.S. Hirata, et al, Phys. Lett. B 205 416 (1988); Phys. Lett. B280 146 (1992). 

[7] 1MB Coll., Casper et al, Phys. Rev. Lett. 66 2561 (1991); Phys. Rev. D46 3720 (1992). 

[8] J. G. Learned and S. Pakvasa, Astropart. Phys. 3, 267 (1995). 

[9] H. Athar, M. Jezabek and O. Yasuda, Phys. Rev. D 62, 10 3007 (2000) [hep-ph/0005104] . 
[10] D. V. Ahluwalia, Mod. Phys. Lett. A 16, 917 (2001) hep-ph/0104316| . 

[11] S. Pakvasa, Mod. Phys. Lett. A 19, 1163 (2004) [Yad. Fiz. 67, 1179 (2004)] [hep-ph/0405179] . 
[12] G. Barenboim and C. Quigg, Phys. Rev. D 67, 073024 (2003) [hep-ph/0301220| . 

[13] J. F. Beacom, N. F. Bell, D. Hoo per, S. Pakvasa and T. J. Weiler, Phys. Rev. D 68, 093005 (2003) [Erratum-ibid. D 72, 

019901 (2005)] [hep-ph/0307025] . 
[14] Y. Chikashige, R. N. Mohapatra and R. D. Peccei, Phys. Rev. Lett. 45, 1926 (1980). G. B. Gelmini and M. Roncadelli, 

Phys. Lett. B 99, 411 (1981); 
[15] A. Acker, S. Pakvasa and J. T. Pantaleone, Phys. Rev. D 45, 1 (1992). 

[16] J. F. Beacom, N. F. Bell, D. Hooper, S. Pakvasa and T. J. Weiler, Phys. Rev. Lett. 90, 181301 (2003) [hep-ph/021 1305 . 
[17] L. Wolfenstein, Nucl. Phys. B 186, 147 (1981); S. T. Petcov, Phys. Lett. B 110, 245 (1982); S. M. Bilenky and B. Pon- 

tecorvo, Sov. J. Nucl. Phys. 38, 248 (1983); C. Giunti, C. W. Kim and U. W. Lee, Phys. Rev. D 46, 3034 (1992); 

M. Kobayashi and C. S. Lim, Phys. Rev. D 64, 013003 (2001). K. R. Balaji, A. Kalliomaki and J. Maalampi, Phys. Lett. 

B 524, 153 (2002). 

[18] J. F. Beacom, N F. Bell, D. Hooper, J. G. Learned, S. Pakvasa and T. J. Weiler, Phys. Rev. Lett. 92, 011101 (2004) 
[hep-ph/030715T] . 

[19] P. Keranen, J. Maalampi, M. Myyrylainen and J. Riittinen, Phys. Lett. B 574, 162 (2003) [hep-ph/0307041| . 



22 



D. Hooper, D. Morgan and E. Winstanley, Phys. Lett. B 609, 206 (2005) [hep-ph/041 0094 . 
M. Gasperini, Phys. Rev. D 39, 3606 (1989). 

H. Minakata and A. Y. Smirnov, Phys. Rev. D 54, 3698 (1996) hep-ph/9601 3"Ti"1; 
R. Fardon, A. E. Nelson and N. Weiner, JCAP 0410, 005 (2004) astro-ph/0309800 . 
P. Q. Hung and H. Pas, Mod. Phys. Lett. A 20, 1209 (2005) |astro-ph/0311131|. 

E. Waxman and J. N. Bahcall, Phys. Rev. Lett. 78, 2292 ( 1997) |astro-ph/970123l] . 

A. Strumia and F. Vissani, Nucl. Phys. B 726, 294 (2005) [hep-ph/0 503246 . 
P. Lipari, Astropart. Phys. 1, 195 (1993). 

W. M. Yao et al. [Particle Data Group], J. Phys. G 33, 1 (2006). 

T. Kashti and E. Waxman, Phys. Rev. Lett. 95, 181101 (2005) [astro-ph/0507599] . 

M. Conversi, E. Pancini and O. Piccioni, Phys. Rev. 71, 209 (1 947). 

J. P. Rachen and P. Meszaros, Phys. Rev. D 58, 123005 (1998) [astro-ph/ 9802280 . 

L. A. Anchordoqui, H. Goldberg, F. Halzen and T. J. Weiler, P hys. Lett. B 593, 42 (2004) [astro-ph/03"T l002 . 
T. Piran, Rev. Mod. Phys. 76, 1143 (2004) [astro-ph/0405503| . 

B. Zhang and P. Meszaros, Int. J. Mod. Phys7A~19, 2385 (2004) ast ro-ph/03113211 . 
A. Dar and A. De Rujula, Phys. Rept. 405, 203 (2004) [astro-ph/0 308248 . 

R.D. Preece et al, Astrophys.J.Suppl.Ser. 126, 19 (2000^"~ 

L. A. Anchordoqui, H. Goldberg, F. Halzen and T. J. Weiler, Phys. Lett. B 621, 18 (2005) [hep-ph/0410003| . 
The details of the Montecarlo code will be presented elsewhere. 
S.L. Glashow, Phys. Rev. 118, 316, (1960). 
A. Roberts, Rev. Mod. Phys. 64, 259 (1992); 
E. Andres et al, Nature 410, 441 (2001). 



A. Achterberg et al. [IceCube Collaboration], Astropart. Phys. 26, 155 (2006) astro-ph/0604450 . 

V. Balkanov et al, Astropart. Phys. 14, 61 (2000). 

S. E. Tzamarias, Nucl. Instrum. Meth. A 502, 150 (2003). 

J. A. Aguilar et al. [ANTARES Collaboration], Astropart. Phys. 26, 314 (2006) [astro -ph/0606229] . 

F. Ameli, M. Bonori and F. Massa, Eur. Phys. J. C 25, 67 (2002). 

J.Learned, Nucl. Phys. Proc. Suppl. 118 (2003) 405. 

I. Kravchenko et al, Int. J. Mod. Phys. A 21S1, 153 (2006). 

H. Falcke, P. Gorham and R. J. Protheroe, New Astron. Rev. 48, 1487 (2004) ast ro-ph/0409229] . 
A. Letessier-Selvon, Nucl. Phys. Proc. Suppl. 118, 399 (2003) [astro-ph/0208526| . 
Euso Colla boration, Nuovo C im. 24C (2001) 445. 

F. Vissani, |astro-ph/0609575| 

D. Meloni and T. Ohlsson, |hep-ph/0612279| 
D. Fargion, |astro-ph/9704205| 

Y. Fukuda et al [Super-Kamiokande Collaboration], Phys. Rev. Lett. 81, 1562 (1998) [hep-ex/9807003] . 
P. D. Serpico an d M. Kachelriess, Phys. Rev. Lett. 94, 211102 (2005) [hep-p h/0502088]. P. Bhattacharjee and N. Gupta, 
|hep-ph/050TT9T| W. Winter, P hys. Rev. D 74, 033015 (2006) |hep-ph/0604191|. P. D. Serp ico, Phys. Rev. D 73, 047301 
(2006) arXiv:hep-ph/0511313 . W. Rodejohann, JCAP 0701, 029 (2007) hep-ph/06 12047 . 



23 




FIG. 1: Distribution of the observable flavor ratio R efl — (ve/v^obs for different models of the v source and different neutrino 
properties. A source model is defined by the relative intensity (summing over v and V) of the emission for the three flavors: 
\y B , Vp, u T ] source- The distribution of _R eM is entirely due to the present uncertainties in the determination of the v mixing 
parameters. The source models [0,1,0], and [1,0,0] correspond to pure or emission; the model [1,1.86,0] corresponds 
to the pion dominated emission from a thin source with a power law spectrum of slope 2. The two dashed curves show the 
distributions of R ell for a j/-decay model |l6j ] where only the lowest v mass eigenstate (in the direct and inverse hierarchy) is 
stable. 




FIG. 2: Distribution of the observable flavor ratio i? MT = (^ M /^ T ) bs (see caption of fig. [T]). Some of the curves have been 
renorrnalized for an easier observation of their shape. 
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FIG. 3: Neutrino flavor ratio (f M +V fl )/(i> < , + z7 e ) at the source obtained from the chain decay of charged pions with a power law 
energy spectrum cx E^ a . Energy losses are assumed to be negligible. The dashed line neglects the effect of muon polarization 
in pion decay. 
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FIG. 4: Neutrino flavor ratio {y^ + u ll ,)/{v E +V e ) at the source for the chain decay of different mesons having a power law 
spectrum cx E~ a . The different lines refer to the decay of it , K and Kl- The dashed (dotdashed) line show the effect of 
neglecting muon polarization in the chain decay of charged pions (kaons). 
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FIG. 5: Characteristic times for a muon plotted as a function of its energy. The different lines correspond to (i) the decay 
time ide C ; (ii) the time tmatter for energy loss in ordinary matter, assuming a gaseous medium with the average composition 
of the interstellar medium (ISM) and the density p = 10~ 10 g cm -3 (this time scales as ^matter p 1 ); (iii) the time t Byn for 
synchrotron losses in a magnetic field B = 10 4 Gauss (t By n oc -B" 1 ), and (iv) the time tic for Inverse Compton losses. The 
target radiation field is assumed isotropic, with energy spectrum of GRB form (|63|l with £b = 1 KeV, and energy density 
p-y — ps (with B = 10 4 Gauss). 




FIG. 6: Characteristic times for charged pions. The different lines correspond to: (i) the decay time idee; (ii) the time tmatter 
for hadronic interactions with a gaseous medium with the ISM composition and density p — 10 -10 g cm" 2 ; (iii) the time t ayn 
for synchrotron losses (in a randomly oriented magnetic field B — 10 4 Gauss); (iv) the time tic for Inverse Compton losses; 
(v) the time t n ^ for photo-hadronic interactions. The last two times are calculated for a radiation field of GRB form (|63[1 with 
Sb = 1 KeV, and energy density p 1 = ps (with B — 10 4 Gauss). Point a (b) is the intersection of the £dcc and t Byn (tmatter) 
curves. The energy that correspond to the point is the critical energy above which synchrotron losses (the hadronic interaction 
probability) are (is) important before decay. For the values of B and p chosen in this illustration the presence of a gas of 
ordinary matter is negligible. 
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FIG. 8: Proton interaction rate in a radiation field of the form (|63(l . The energy E* (defined in equation l|72p) corresponds 
to the threshold proton energy for interaction with photons with energy Sh- The dashed line neglects the interactions with all 
photons with energy e > Eh- 
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FIG. 9: Distribution of e r (the photon energy in the proton rest frame) for the p7 interactions of protons of different energy 
E p . The isotropic photon field has the energy spectrum of GRB form (I63|) . The dashed line neglect the contribution of photons 
above the break energy Eh- 




FIG. 10: Energy distributions (in the form Qj(E)E 2 versus E) for different secondary particles (itr^, K ± , Kl) produced by 
the P7 interactions of relativistic protons in a GRB radiation field. The shape of the primary proton spectrum is also shown 
(in the form N P (E) E 2 versus E) as a thin solid line. The energy spectrum of the isotropic target photons has the form (|63p . 
All energies are in units of the characteristic energy E* (defined in equation J72}). The dashed lines are the distributions of 
the meson energy at decay assuming the presence of magnetic field with value B = 2.78 x 10 3 (eb/KeV) Gauss. 
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FIG. 11: Neutrino energy spectrum (summed over all flavors) from a GRB source, according to the model [25|]. The calculation 
is performed assuming a proton spectrum N P (E) cx E~ a , traveling in an isotropic radiation field of form (I63|l , The different 
curves are obtained varying the exponent a of the proton spectrum, or the exponents /3i and $2 of the radiation field. All 
energies are measured as fraction of E* = m p eth/(2eb). The energy losses of all secondary particles are assumed to be negligible. 
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FIG. 12: Neutrino energy spectrum (summed over all flavors) from a GRB source (see text). The primary proton spectrum 
has slope a — 2, while the target photons have a broken power law energy distribution (see equation (JB3J)) with slopes /3i = 1 
and P2 = 2 below and above the break energy e\>. The curves are calculated assuming different values for the magnetic field, 
corresponding to values of the critical energy for muons e M = 00, 30 and 3 (see eq. ((74)) '). 
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FIG. 13: Energy spectra for and F e from a GRB source (proton slope a = 2, target photon slopes /3i = 1 and fa = 2, 
no magnetic field) For the F e flux also the contributions from meson and neutron decay are separately shown. Neutron decay 
dominates at low energy. 
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FIG. 14: Ratio (u^ + v^)/(y e + v<-) for neutrinos from a GRB source. The ratio correspond to the spectra shown in fig. 1111 
The rapid increase at low energy is due to the effect of 77 e from n decay. 
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FIG. 15: Ratio (y^ + V tl )/{y e +V e ) for neutrinos from a GRB source. The ratio correspond to the spectra shown in fig. 1121 
The three curves are calculated assuming different values of the magnetic field in the source (and proton slope a — 2, target 
photon slopes j3i = 1 and @2 = 2). The decrease of the flavor ratio at high energy for large B is associated to the important 
synchrotron losses of high energy muons. 
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FIG. 16: Ratio (v e /V e ) at the source for neutrinos from a GRB source. The ratio correspond to the spectra shown in fig. 1121 
At low energy the ratio is low because of the importance of V e from n decay. At higher energy the ratio decreases because of 
the increased importance of 7r~ production. For high magnetic field, at high energy the dominant source of u e and V e becomes 
Kl decay, and the (v e /V e ) ratio becomes unity. 
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FIG. 17: Ratio (f M /F M ) at the source for neutrinos from a GRB source. The ratio correspond to the spectra shown in fig. 1121 
The three curves correspond to three different values of the magnetic field in the source. At low energy (and at all energies for 
small B) the ratio is unity, reflecting the fact that the chain decay of a charged pion produces a and a 77 M . In the presence of 
a large magnetic field, the neutrinos from muon decay are suppressed, and the flavor ratio is determined by the ir + /n~ ratio. 
At the highest energies charged kaon decay becomes the dominant source of and and the flavor ratio reflects the K + / K~ 
ratio. 




FIG. 18: Observable ratio (yn + V ll )/(y e + Ve) after propagation for the neutrinos from a GRB source. The solid lines 
corresponds to the source models shown in fig. 1121 and 1151 The three curves differ for the value of the magnetic field that 
corresponds to e M = 3, 30 and oo (no field). The neutrino oscillation parameters have been chosen as #23 = 45°, #12 = 34° and 
#13 = 0. The thin dotted are calculate for the "naive" assumption of a source composition [i/ E ,v l _ l ,i> T ] = [1,2,0]. The dotted 
line assumes the same mixing parameters as above. The Dashes (dot-dashed) line has 823 = 50° ( 40°). 



